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We have investigated atom-nanotube van der Waals (vdW) coupling in atomically doped carbon 
nanotubes (CNs). Our approach is based on the perturbation theory for degenerated atomic levels, 
thus accounting for both weak and strong atom- vacuum-field coupling. The vdW energy is described 
by an integral equation represented in terms of the local photonic density of states (DOS). By solving 
it numerically, we demonstrate the inapplicability of standard weak-coupling-based vdW interaction 
models in a close vicinity of the CN surface where the local photonic DOS effectively increases, giving 
rise to an atom-field coupling enhancement. An inside encapsulation of atoms into the CN has been 
shown to be energetically more favorable than their outside adsorption by the CN surface. If the 
atom is fixed outside the CN, the modulus of the vdW energy increases with the CN radius provided 
that the weak atom-field coupling regime is realized (i.e., far enough from the CN). For inside atomic 
position, the modulus of the vdW energy decreases with the CN radius, representing a general effect 
of the effective interaction area reduction with lowering the CN curvature. 

PACS numbers: 61.46.+W, 73.22.-f, 68.65.-k, 34.50.Dy 
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I. INTRODUCTION 

Carbon nanotubes (CNs) are graphene sheets rollcd- 
up into cylinders of approximately one nanometer in di- 
ameter. Extensive work carried out worldwide in recent 
years has revealed the intriguing physical properties of 
these novel molecular scale wires [J, Q ■ Nanotubes have 
been shown to be useful for miniaturized electronic, me- 
chanical, electromechanical, chemical and scanning probe 
devices and materials for macroscopic composites [3| . Im- 
portant is that their intrinsic properties may be substan- 
tially modified in a controllable way by doping with ex- 
trinsic impurity atoms, molecules and compounds Q|. Re- 
cent successful experiments on encapsulation of single 
atoms into single-wall CNs |5| and their intercalation 
into single- wall CN bundles UH, along with numerous 
studies of monoatomic gas absorbtion by the CN bun- 
dles (see for a review) , stimulate an in-depth analysis 
of the atom-CN van der Waals (vdW) interactions. The 
problem is both of fundamental and of applied interest, 
because it sheds light on the peculiarities of the atom- 
electromagnetic-field interactions in low-dimensional dis- 
persive and absorbing surroundings, on the one hand, 
and, on the other, the atom-nanotube vdW interaction 
is the most important one determining the filling factor 
in encapsulating single atoms into CNs. 

The vdW interaction is the result of electrodynamic 
coupling between the polarization states of the interact- 
ing entities that results from the vacuum fluctuations 
of the electroma gnet ic field. A theoretical analysis re- 
cently done H, ffl. Il0| of spontaneous decay dynamics of 
excited atomic states near CNs has brought out fasci- 
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nating peculiarities of the vacuum-field interactions in 
atomically doped CNs. In particular, close to the nan- 
otube, the relative density of photonic states (DOS) effec- 
tively increases due to the presence of additional surface 
photonic states coupled with CN electronic quasiparticle 
excitations. This causes an atom-vacuum-field coupling 
constant, which is proportional to the photonic DOS, to 
be very sensitive to an atom-CN-surface distance. If the 
atom is close enough to the CN surface and the atomic 
transition frequency is in the vicinity of a resonance of the 
photonic DOS, the system shows strong atom-vacuum- 
field coupling giving rise to rearrangement (" dressing" ) of 
atomic levels by the vacuum-field interaction. If the atom 
moves away from the CN surface, the atom-field coupling 
strength decreases, smoothly approaching the weak cou- 
pling regime at large atom-surface distances since the 
role of surface photonic states diminishes, causing the 
relative photonic DOS to decrease. This suggests strictly 
nonlinear atom-field coupling and a primary role of the 
distance-dependent (local) photonic DOS in the vicin- 
ity of the CN, so that (weak-coupling) vdW interaction 
models based upon standard vacuum quantum electrody- 
namics (QED) as well as those using the linear response 
theory (see [Tllll2| for a review) are, in general, inappli- 
cable for an atom in a close vicinity of a carbon nanotube. 

To give this issue a proper quantum electrodynamic 
consideration and thereby to advance in physical under- 
standing the vdW interactions deeper than simple ana- 
lyzing empirical model potentials [13|, we develop a uni- 
versal quantum approach to the vdW energy calculation 
of a (two-level) atomic system in the vicinity of an in- 
finitely long single-wall CN. The approach is based upon 
the perturbation theory for degenerate atomic levels (see, 
e.g., [13 ) j thereby allowing one to account for both strong 
and weak atom-vacuum-field coupling regimes. We derive 
an integral equation for the vdW energy of the two-level 
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atomic system near the CN. The equation is represented 
in terms of the local photonic DOS and is valid for both 
strong and weak atom-field coupling. By solving it nu- 
merically, we demonstrate the inapplicability of weak- 
coupling-based vdW interaction models in a close vicin- 
ity of the nanotube surface. We also show that an inside 
encapsulation of doped atoms into the nanotube is ener- 
getically more favorable than their outside adsorption by 
the nanotube surface. 

The rest of the paper is arranged as follows. Section ITT1 
presents a theoretical model we use to derive a univer- 
sal integral equation for the ground-state vdW energy of 
the two-level atomic system interacting with a vacuum 
electromagnetic field in the vicinity (inside or outside) 
of the CN. In describing the atom-field interaction, we 
follow an original line of Ref. [TEl flol ]. adopting their 
electromagnetic field quantization scheme for a particu- 
lar case of the field near an infinitely long achiral single- 
wall CN. The carbon nanotube is considered to be an in- 
finitely thin anisotropically conducting cylinder. Its sur- 
face conductivity is represented in terms of the 7r-electron 
dispersion law obtained in the tight-binding approxima- 
tion with allowance made for azimuthal electron momen- 
tum quan tization and axial electron momentum relax- 
ation 13 ■ Only the axial conductivity is taken into ac- 
count and the azimuthal one, bein g st ron gly supp ressed 
by transverse depolarization fields |l8t Il9l l2Ct \2l\ . is ne- 
glected. In Sec. IIIII the vdW energy equation derived in 
Sec. |n] is analyzed qualitatively to demonstrate its va- 
lidity and correctness in both weak and strong atom- 
field coupling regimes. We show that the equation ob- 
tained reproduces a well-known perturbation theory re- 
sult in the weak coupling regime and reduces to the inte- 
gral equation of a simplified form in the strong coupling 
regime. Section Hvl presents and discusses the results of 
the various numerical calculations of the vdW energies of 
the two-level atom placed inside or outside achiral CNs 
of different radii. A summary and conclusions of the work 
are given in Sec. Ivl 



II. THE GROUND-STATE VAN DER WAALS 
ENERGY 

A. The field quantization formalism 

The quantum theory of the van der Waals interaction 
between an atomic system and a carbon nanotube in- 
volves an electromagnetic field quantization procedure in 
the presence of dispersing and absorbing bodies. Such a 
procedure faces difficulties similar to those in quantum 
optics of 3D Kramers-Kronig dielectric media where the 
canonical quantization scheme commonly used does not 
work since, because of absorption, correspondin g op era- 
tor Maxwell equations become non-Hermitian [22| . As 
a consequence, their solutions cannot be expanded in 
power orthogonal modes and the concept of modes it- 
self becomes more subtle. We, therefore, use a unified 




FIG. 1: (Color online) The geometry of the problem. 



macroscopic QED approach developed in Refs. \pl , Il6| 
and adapted for an atom near a CN in Refs. pi llfij. In 
this approach, the Fourier-images of electric and mag- 
netic fields are considered as quantum mechanical observ- 
ables of corresponding electric and magnetic field opera- 
tors. The latter ones satisfy the Fourier-domain operator 
Maxwell equations modified by the presence of a so-called 
operator noise current density J(r, to) written in terms of 
a three-dimensional (3D) vector bosonic field operator 
i(r,uj) and a medium dielectric tensor e(r,ui) (supposed 
to be diagonal) as 

ii( r > w ) = — v / ftIm^i(iV^) /i( r »^)) i = 1,2,3. (1) 

Z7T 

This operator is responsible for correct commutation re- 
lations of the electric and magnetic field operators in the 
presence of medium-induced absorbtion. In this formal- 
ism, the electric and magnetic field operators are ex- 
pressed in terms of a continuum set of the 3D vector 
bosonic fields f(r, oj) by means of the convolution over r 
of the current (JJJ with the classical electromagnetic field 
Green tensor of the system. The bosonic field operators 
P(r,oj) and f(r, oj) create and annihilate single-quantum 
electromagnetic medium excitations. They are defined by 
their commutation relations 

[/i(r,w),/J(rV)] = M(r-r')^-^'), (2) 

[fau), = [#(r, W ),/7(r>')] =0 

and play the role of the fundamental dynamical variables 
in terms of which the Hamiltonian of the composed sys- 
tem " electromagnetic field + dissipative medium" is writ- 
ten in a standard secondly quantized form as 

H = Jduhufdrit(T,u)-f(T,w). (3) 

Consider a neutral atomic system with its centre-of- 
mass positioned at the point ya near an infinitely long 
single-wall CN. Since the problem has a cylindric sym- 
metry, it is convenient to assign the orthonormal cylin- 
dric basis {e r , e v , e z } in such a way that e z is directed 
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along the nanotube axis and ta — Va^t — {ta, 0, 0} (see 
Fig. nj. For carbon nanotubes, their strong transverse 
depolarization along with transverse electron momentum 
quantization allow one to neglect the azimuthal current 
and radial polarizability |l7lll8lll9ll2Ctl21j . in which case 
the dielectric tensor components e rr and e vv are identi- 
cally equal to unit. The component e zz is caused by the 
CN longitudinal polarizability and is responsible for the 
axial surface current parallel to the e z vector. This cur- 
rent may be represented in terms of the ID bosonic field 
operators by analogy with Eq. (JIJ. Indeed, taking into 
account the dimensionality conservation in passing from 
bulk to a monolayer in Eq. and using a simple Drudc 
relation [Hi 



s (R,w) 



e zz (R,m) - 1 



(4) 



where R = {R cn , 0, Z} is the radius- vector of an arbitrary 
point of the CN surface, R cn is the radius of the CN, 
a zz (R, to) is the CN surface axial conductivity per unit 
length, S is the area of a single nanotube, px is the tubule 
density in a bundle, one immediately has from Eq. Q 



t,d , hujRea zz (R,u!) - 
J(R, ui) = \ /(R, oj)e z 



(5) 



with /(R, oj) being the ID bosonic field operator defined 
on the CN surface. The total nonrelativistic Hamilto- 
nian of an atomic system interacting with a CN modi- 
fied quantized vacuum electromagnetic field can then be 
represented in the form 



H= / diuhuj dRf(R,Lu)f(R,uj) 



E — 



Ivrii 



Pi -A(fj) 

c 



+ 2 / dr PA(r)<p A {r) + I dr p a (t)<p(t) . 



(6) 



where wii, qi, fj and pi are, respectively, the masses, 
charges, coordinates (relative to r^) and momenta of 
the particles constituting the atomic subsystem. The 
first term of the Hamiltonian describes the medium- 
assisted vacuum electromagnetic field in the presence of 
the CN. The second and the third terms of the Hamilto- 
nian represent the kinetic energy of the charged particles 
and their mutual Coulomb interaction, respectively, with 



being the scalar potential and 

PA(r) =J2q l 5(r-r i 



(7) 



(8) 



being the charged density of the atomic subsystem. The 
last term accounts for the Coulomb interaction of the par- 
ticles with the CN. The vector potential A and the scalar 



potential (p of the CN-modified electromagnetic field are 
given for an arbitrary r = {r, tp, z} in the Schrodingcr 
picture by 



(9) 



A(r) = / dio—'E (r,w)+h.c, 



-V<p(r)= / dwE (r,w)+h.c, 



(10) 



where 



E ±(ll) (r,w) =/'dr / * J -^(r-r')-i(r / ,w) (11) 

is the transverse (longitudinal) electric field operator 
with 



and 



1 

47rr 



5^(r)=6 af3 8(r)-6Ur) 



(12) 



(13) 



being the longitudinal and transverse dyadic (5-functions, 
respectively, and E representing the total electric field op- 
erator which satisfies the following set of Fourier-domain 
Maxwell equations 



V x E(r,w) = ijfeH(r,£j), 



47T- 



(14) 



V x H(r, u>) = -ifeEfr, w) H I(r, cA (15) 

c 

Here, H stands for the magnetic field operator, k = uj/c, 
and 



1(t,lu) = JdR5(r~R) J(R,w) 

= 2J(R cn ,ip,z,u)6(r-R cn ), (16) 

is the exterior operator current density [with J(R, u>) de- 
fined by Eq. JSJl] associated with the presence of the CN. 
From Eqs. in"4l- ((Pojl in view of Eq. JSJ, it follows that 

4tt f 

E(r,w) =%— k /dRG(r,R,w)-J(R,w) (17) 



[and H = (i/c) _1 V x E accordingly], where G is the 
Green tensor of the classical electromagnetic field in the 
vicinity of the CN. The set of Eqs. (JBJ-JTJJl forms a closed 
electromagnetic field quantization formalism in the pres- 
ence of dispersing and absorbing media which meets all 
the basic requirements of a standard quantum electro- 
dynamics [16J . All information about medium (a CN in 
our case) is contained in the Green tensor G via the CN 
surface axial conductivity in Eq. (J5J. The Green tensor 
components satisfy the equation 

(VxVx-^G M (r,R, w )={(r-R), (18) 

a—r,ip,z 
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together with the radiation conditions at infinity and the 
boundary conditions on the CN surface. This tensor was 
derived and analysed in Ref. [Tol | . 

Assuming further that the atomic subsystem is suffi- 
ciently localized in space, so that the long-wavelength 
approximation applies, one can expand the field oper- 
ators A(r) and <p(r) in the Hamiltonian © around the 
atomic center of mass position ya and only keep the lead- 
ing non-vanishing terms of the expansions. Then, under 
the condition of the Coulomb gauge [pi, A] — 0, one ar- 
rives at the approximate total Hamiltonian of the form 



H = H F + H A + H 



(i) 

AF 



H 



(2) 
AF i 



(19) 



where 



H F = 



w) , (20) 



(21) 



H 



(i) 

AF 



ITliC 



p r A(r A ) + d-V<p(r A ) , 



H {2) - 

AF — 



E 



2rrii c 



' A 2 (r A ) 



(22) 



(23) 



are, respectively, the Hamiltonian of the vacuum electro- 
magnetic field modified by the presence of the CN, the 
Hamiltonian of the atomic subsystem, and the electric 
dipole approximation for the Hamiltonian of the atom- 
field interaction (separated into two contributions ac- 
cording to their role in the vdW energy - see below) with 
d = qiTi being the electric dipole moment operator 
of the atomic subsystem. 



B. The two- level approximation 

Starting from Eqs. JSJ-lOJl and using Eqs. ©-(O, 
(|17|l and (JSJ , one obtains the following two- level approx- 
imation for the total Hamiltonian of the system (see Ap- 
pendix 



H = Tio + Hi 



(24) 



where the first term stands for the unperturbed Hamil- 
tonian given by 



Ho = H F + Ha 



(25) 



with Hp being the Hamiltonian (|20|l of the CN modified 
field and 



Ha = —r~ c 2 



(26) 



representing the 'effective' unperturbed two-level atomic 
subsystem, and the second term stands for their interac- 
tion given by 

H int =y < Ly'dR[gW(r A) R,a;)at 

-- g M (r A ,R, W )a]/(R, W ) + h.c. (27) 
Here, the Pauli operators 

<r=\l)(u\, tf = \u){l\, o z = \u)(u\-\l)(l\, (28) 
\u){u\ + \l){l\ = I, 

describe electric dipole transitions between the two ato- 
mic states, upper \u) and lower separated by the tran- 
sition frequency u>a- This 'bare' transition frequency is 
modified by the interaction which, being indepen- 
dent of the atomic dipole moment, does not contribute 
to mixing the \u) and |/) states, giving rise, however, to 
the new renormalized transition frequency 



UJA = OJA 



(TlU> A )' 



^dujJdR\ g x (T A ,n,uj) 



(29) 



in the 'effective' atomic Hamiltonian (|26|l . On the con- 
trary, the interaction l|22l) . being dipole moment depen- 
dent, mixes the \u) and |^) states, yielding the perturba- 
tion Hamiltonian l|27|l with the (dipole) interaction ma- 
trix elements of the form 

gW (r A , R, w) = g X (r A , R, «) ± — g" (r a , R, w) , (30) 

UJ A 



d z ^nhuj Re o zz (R, uj) ±{ ^G zz {r A , R, w) (31) 



where 

Au>A 

with 

^G zz (r A ,-R,uj) = fdrSj-PitA-T) G zz (r,K,uj) 

(32) 

being the field Green tensor zz-component transverse 
(longitudinal) with respect to the first variable. This is 
the only Green tensor component we have to take account 
of. All the other components can be safely be neglected, 
because our model neglects the azimuthal current and 
radial polarizability of the CN from the very outset. 

The matrix elements (|30[) and (|31[) have the following 
properties (see Appendix |B|) 



JdR\g ±{in (r A ,-R,uj) 



(Huja) 



2ttuj 



(33) 



and 



/'dR|g( ± )(r A ,R,«)| a = 



UJ A 



{h " A)2 r (uj) 



Z ll (rA,u) 



(34) 
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Here, £ x(| D(r,i , oS) is the transverse (longitudinal) local 
photonic DOS function defined by |10l| 



where 



T^{r A ,u>) = ^dllm^G^){v M v A ,u>) (36) 

is the transverse (longitudinal) atomic spontaneous decay 
rate near the CN with 

^GjrP(r A ,r A ,u) = J drdr'5^\r A - r) 

x G zz (r,r', W )^ (ll) (r'-r A ) (37) 

being the Green tensor zz-component transverse (longi- 
tudinal) with respect to both variables, and Tq(uj) is the 
same quantity in vacuum where |23| 



(35) 



ImG 



6irc 



(38) 



Note that the transverse and longitudinal imaginary 
Green tensor components in Eq. (|36|l can be written as 

Im X G x = ImGL+Im^, (39) 
LnllGf, = ImllGl (40) 

with - L 0I)G^"^ representing the "pure" CN contribu- 
tion to the total imaginary Green tensor. The longitu- 
dinal imaginary Green tensor in Eq. (|40|l is totally con- 
tributed by a (longitudinal) CN static polarization field 
and, therefore, does not contain the vacuum term which 
is the transverse one by its definition as there are no po- 
larization Coulomb sources in vacuum. Thus, Eq. Ij35(l 
may, in view of Eqs. Ij36|l - (|40|l . be rewritten in the form 



^(t a ,lj) = l + £ (r A ,w) 
£"(rA,w) = £ (r A ,u), 



(41) 
(42) 



where the r^-dependent terms come totally from the 
presence of the CN. 

Eqs. H24|l - 131(l represent the total secondly quantized 
Hamiltonian of an atomic system interacting with the 
electromagnetic field in the vicinity of the carbon nan- 
otube. In deriving this Hamiltonian, there were only 
two standard approximations done. They are the elec- 
tric dipole approximation and the two-level approxima- 
tion. The rotating wave approximation commonly used 
was not applied, and the (diamagnetic) A 2 -term of the 
atom-field interaction [Eq. J23J] was not neglected. The 
latter two approximations are justified in analyzing an 
atomic state evolution due to real dipole transitions to 
the nearest states [24| , described for two-level systems by 
the g( + )-terms of the Hamiltonian l|27|l . They are, how- 
ever, inappropriate in calculating the energy, where the 
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FIG. 2: (Color online) Dimensionless unperturbed energy lev- 
els [given by the Hamiltonian 125H 1 and total ground-state 
energy [given by Eqs. I|46^ - l|50^ ] as functions of the atomic 
position for the two-level atom outside the (9,0) CN. 



initial energy level is being shifted due to the perturba- 
tion caused by virtual dipole transitions via intermediate 
atomic states, so that the system always remains in the 
initial state and the terms neglected contribute to the 
level shift to the same order of magnitude. For two-level 
systems, this is described by g(~'-terms of the Hamil- 
tonian H27|) (usually neglected within the rotating wave 
approximation) together with the atomic transition fre- 
quency l|29ll renormalized due to the (usually discarded) 
A 2 -term of the atom-field interaction. 



C. The van der Waals energy 

The Hamiltonian l|24|) - (|27|l is a starting point for the 
atom-nanotube vdW energy calculation. Supposing the 
two-level atom to be in the ground state, the vdW en- 
ergy is nothing but the r^-dependent contribution to the 
ground-state eigenvalue of this Hamiltonian. Important is 
that the unperturbed atomic subsystem is now described 
by the Hamiltonian (12(it with the renormalized transition 
frequency l|29|l which is represented in terms of the local 
photonic DOS as 



LO A = LO A 



1 

71" ./, 



dw — ^^(r A ,uj) 



(43) 



As is seen from Eq. Q43fl , the transition frequency Q A de- 
creases with increasing the transverse local photonic DOS 
^(ta.ui), i.e. when the atom approaches the CN sur- 
face [lOj , thereby bringing the unperturbed atomic levels 
together, or even making them degenerated if £ (r A) lo) 
is large enough (see a typical example in Fig. |2J. If one 
uses the standard perturbation theory, this would yield 
the divergence of high-order corrections to the vdW en- 
ergy at small atom-CN-surface distances because of the 
small energy denominators in the corresponding power- 
series expansions. To account for a possible degeneracy of 
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the unperturbed atomic levels in a correct way, one has 
to use the perturbation theory for degenerated atomic 
levels (see, e.g., 0]). In so doing, one has to take the 
upper state degeneracy into account of the whole system 
with respect to R, so that the ground-state wave function 
of the whole system should be represented as a coherent 
mixture of the lower and upper states of the form 

|V) = C l \l)\{0}) + / L/dRC u (R,w)|u)|{l(R,w)}>, 



(44) 

where |{0}) is the vacuum state of the field subsystem, 
|{1(R, w)}) is its single-quantum excited state, C\ and 
C U (R, lu) are unknown mixing coefficients of the lower 
and upper states of the whole system. 

The total energy E of the ground state is now given by 
the solution of a secular equation obtained by applying 
the Hamiltonian J2U)-(I23 to the wave function This 
yields the integral equation 



E 



Tiuja 



duj / cZR 



fiU) ■ 



Tiuja 



(45) 



- E 



which the ground-state-atom vdW energy E vw {ya) is de- 
termined from by means of an obvious relation 



E 



Tiuja 



+ E vw (r A ) ■ 



(46) 



with the first term representing the 'bare' (non-inte- 
racting) two-level atom in the ground-state. Using di- 
mensionless variables 



E vw (y a ) 
270 



fir (x) 

270 



7T- (47) 
27o 



with 70 = 2.7 eV being the carbon nearest neighbor hop- 
ping integral appearing in the CN surface axial conduc- 
tivity a zz in Eq. JHJ|, and Eqs. {2}, l|H )l -l(4*3 | . one can 
represent the dimensionless vdW energy e vw (ya) in terms 
of the transverse and longitudinal distance-dependent 
photonic DOS's in the vicinity of the nanotube. In so 
doing, one arrives at the integral equation of the form 



1 \ _ A 

ZlT 



[•OC 



Jo x 



2tt 



<ir 



fo(a:) 



£ (r A ,x) + ( — ) £ (r A ,x) 



x + xa 



(48) 



with £ ^ (*a, x) given for ta > R C n by ( Appendix lB|l 



£ (*A,X) 



3 ^ r c 



dy s (R cn ,x)v(ypl 2 [v {y)u(R cn ) x] K 2 [v {y)u{r A )x\ 
s(R cn , x)v(y) 2 Ip[v(y)u(R cn )x\Kp[v(y)u(R cn )x] 



(49) 



C"(r A ,a;) 



3 ^ [di 



dys(R cn ,x)y 2 v(y) 2 I 2 [v{y)u(R cn )x]K 2 [v{y)u{r A ) 



P— — OO 



s(R C n, x)v(y) 2 I p [v(y)u(R cn )x]K p [v(y)u(R cn )x] 



(50) 



where I p and K p are the modified cylindric Bessel func- 
tions, v(y) = yjy 2 - 1 , u(r) = 2~/ r/Tic, and s{R cn ,x) = 
2iau(R cn )xa zz (R cn ,x) with a zz = 2irh<7 zz /e 2 being the 
dimensionless CN surface conductivity per unit length 
and a — e 2 /hc — 1/137 representing the fine-structure 
constant. The integration contour C goes along the real 
axis of the complex plane and envelopes the branch points 
y = ±l of the function v(y) in the integrands from below 
and from above, respectively. For ta < Rem Eqs. (|49|l and 
(|50|l are modified by a simple replacement ta R C n in 
the Bessel function arguments in the numerators of the 
integrands. 



III. QUALITATIVE ANALYSIS 



Eqs. (|48() - (|50|l describe the ground-state-atom vdW 
energy near an infinitely long single-wall carbon nan- 
otube in terms of the local (distance-dependent) photonic 
DOS. Eq. (|48|l is universal in the sense that it covers both 
strong and weak atom-field coupling regimes defined as 
those violating and non- violating, respectively, the appli- 
cability domain of the conventional stationary perturba- 
tion theory |40j. 
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Weak coupling regime 



expect the condition 



If, when the atom is faraway from the CN, the local 



photonic DOS £ (r^, x) is such small that 



2tt 



and 



°, To(x) 7 i, , 
dx-^C (rx,ar)<l 



(51) 



(52) 



then Eq. I|48|l yields a well-known perturbation theory 
result (see, e.g., .11]) of the form 



XA_ 

2tt J 



dx 



fo(aO 7 



£ (rA,a?) 



2vr 



dx 



~ £ (r A ,x) + — K (rA,a;) 

r n (x) \xaJ 

X 2 X + XA 



, (53) 



where the first term comes from the unperturbed Hamil- 
tonian l|25[l . i|26|) and the second one is the second or- 
der correction due to the perturbation (|27|l . This can be 
equivalently rewritten in the form 



27T J x(x + XA) 



X -p\\, , 
XA 



(54) 



from which, in view of Eqs. (|49|l . H5Q|I and basic properties 
of the modified cylindric Bessel functions [26j, one can 
immediately come to an interesting conclusion. Namely, 
if the atom is fixed outside (inside) the CN in a way 
that the weak atom-field coupling regime is realized, 
i.e. far enough from the CN surface, then the modulus 
of the atom-nanotube vdW energy increases (decreases) 
with the CN radius. The conclusion is physically clear 
as the effective atom-nanotube interaction area is larger 
(smaller) for larger-radius nanotubes when the atom is 
outside (inside) the CN. Important, however, is that this 
obvious conclusion is, strictly speaking, only valid in the 
weak atom-field coupling regime for the outside atomic 
position, while for the inside atomic position it represents 
a general effect of the effective interaction area reduction 
with lowering the CN surface curvature. 

In the large CN radius limit, Eq. (|54|l along with 
Eqs. (|49|) and (|50|l can be shown to reproduce a well- 
known Casimir-Polder result |27] ] for an atom near an 
infinitely conducting plane (see Appendix ICf . 



2. Strong coupling regime 

When the atom is close enough to the CN surface, the 
local photonic DOS £ (r A ,x) is large, so that one might 



1 

2~7T Jo 



°°dx T -^t(rA,x)~l 



(55) 



to hold true. Under this condition, Eq. (|48|l reduces to 
an integral equation of the form 



2n 



dx 



n 



£ {*A,X) 



r 2 r°° 

C _A I 

2n Jo 



£ (r A ,x) + ( — J £ (r A ,x) 

XA 



x - e vw {vA) 



(56) 



valid in the strong atom-field coupling regime. 

Eq. (|55|l corresponds to the inverted unperturbed 
atomic levels described by the Hamiltonian (|26|l with uja 
given by Eq. I|43(l (see Fig. as an example). The levels 
are degenerated when the left-hand side of Eq. is 
exactly 0.5. A noteworthy point is that if one used the 
weak coupling approximation of Eq. I|48l) for the atom 
close to the CN surface where Eq. I)51|l is not satisfied 
and Eq. (|55l) holds true instead, then the result would 
be divergent at the lower limit of integration. To avoid 
this fact, one has to use either strong coupling approxi- 
mation (|56|l or the exact equation l|48|l to calculate the 
vdW energy of the atom close to the nanotube surface. 



IV. NUMERICAL RESULTS AND DISCUSSION 

Using Eqs. I|46 fl -I|5l) |) . we have simulated the total 
ground state energy of the whole system "two-level 
atom + CN-modified vacuum electromagnetic field" and 
the vdW energy of the two-level atom nearby metal- 
lic and semiconducting carbon nanotubes of different 
radii. The dimensionless free-space spontaneous decay 
rate in Eq. (|48|l was approximated by the expression 
To (a;) w Aa 3 x 3 /3x 2 i (a = 1/137 is the fine-structure con- 
stant) valid for atomic systems with Coulomb interac- 
tion Q . 

Figure 21 shows the dimensionless ground state energy 
level of the whole system [given by the corresponding 
eigenvalue of the the total Hamiltonian l|19|) and repre- 
sented by Eqs. 1)46 - l|5L) p ] and two dimensionless energy 
levels of the unperturbed Hamiltonian l|25l) , 126(1 as func- 
tions of the atomic position outside the (9,0) CN. As the 
atom approaches the CN surface, its unperturbed levels 
come together, then get degenerated and even inverted 
at a very small atom-surface distance. In so doing, the 
weak coupling approximation for the ground state energy 
[given by Eqs. gHJ, 63, and g5J, J5UJ] diverges 
near the surface, whereas the strong coupling approxima- 
tion [Eqs. gBJl, (03), (2| and gSJ, JSUJl] yields a finite 
result. The exact solution reproduces the weak-coupling 
approximation at large and the strong-coupling approxi- 
mation at short atom-surface distances, respectively. 
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(a) (9,0) 
r = 1.5R ; r-2.5H 

A at ' 

— - r =3.5R 



1w 




0.330.520.58 




exact solution 

strong coupling approx. 

weak coupling approx. 



FIG. 3: (Color online) Transverse local photonic DOS's (a) 
and ground-state vdW energies (b) as functions of the atomic 
position for the two-level atom outside the (9,0) CN. The 
'bare' atomic transition frequencies are indicated by dashed 
lines in Fig. EJa). 



The degeneracy of the unperturbed atomic levels is 
the consequence of the transverse local photonic DOS 
increase [see Eq. as the atom approaches the CN 

surface. A typical example is shown in Fig. (a) for 
the atom outside the (9,0) CN. The transverse DOS 
function is seen to increase with decreasing the atom- 
CN-surface distance, representing the atom-field coupling 
strength enhancement with surface photonic modes of 
the nanotube as the atom approaches the nanotube sur- 
face |g, [lOj . The vertical dashed lines indicate the 'bare' 
atomic transition frequencies xa for which the vdW en- 
ergies shown in Fig. 0(b) were calculated. Some of them 
are the positions of peaks, the others are the positions 
of dips of the local photonic DOS. They are typical for 
some rear-earth ions such as Europium, or for heavy 
hydrogen-like atoms such as Caesium (supposed to be 
non-ionized near the CN). For example [28| . for a well- 
known two-level dipole transition 7 F% <-> 5 Do of Eu- 



ropium at A = 615 nm one has xa ~ 0.37, whereas 
for Caesium one obtains from its first ionization poten- 
tial the estimate x A » 3.89 eV x 3/4 x (270)" 1 » 0.5 
(the factor 3/4 comes from the Lyman series of Hydro- 
gen), or less for highly excited Rydberg states. 

Figure |3(b) shows the vdW energies of the atom out- 
side the (9,0) CN for three different 'bare' atomic tran- 
sition frequencies [indicated by vertical dashed lines in 
Fig. El (a)]. Here, exact numerical solutions of Eqs. 1(35)1- 
(|50|l are compared with those obtained within the weak 
coupling approximation ((54(1 ■ l|49(l . 1 (50(1 and within the 
strong coupling approximation 1(56(1 . 1(49(1 . 1(50(1 . At small 
atom-CN-surface distances, the exact solutions are seen 
to be fairly well reproduced by those obtained in the 
strong coupling approximation, clearly indicating the 
strong atom-field coupling regime in a close vicinity of the 
nanotube surface. The deviation from the strong coupling 
approximation increases with transition frequency xa. 



(a) (10,10) 

r = 0; r=G2R a ; 

r = 0.6R 




0.220.30.45 1 



r./R 

A cn 



0,0 0,2 0,4 0,6 1,4 
0,0 -r— * ' ' 1 ■ *—// — L - 



1,6 l,g 



04 




(b) (10,10) 
exact solution 
x =0.22; 

A 9 

x=0.30; 

•• x=0.45 



FIG. 4: (Color online) Longitudinal local photonic DOS's (a) 
and ground-state vdW energies (b) as functions of the atomic 
position for the two-level atom nearby the (10,10) CN. The 
'bare' atomic transition frequencies are indicated by dashed 
lines in Fig. a). 
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that is easily explicable since the degeneracy condition 
(|55|1 of the unperturbed atomic levels is more difficult 
to reach for larger inter-level separations. As the atom 
moves away from the CN surface, the exact solutions de- 
viate from the strong coupling solutions and approach 
those given by the weak coupling approximation, indicat- 
ing the reduction of the atom-field coupling strength with 
raising the atom-surface distance. The weak-coupling so- 
lutions are seen to be divergent close to the nanotube 
surface as it should be because of the degeneracy of the 
unperturbed atomic levels in this region. 

Figure (a) shows the longitudinal local photonic 
DOS's for the atom inside the (10,10) CN at three repre- 
sentative distances from the nanotube wall. As the atom 
approaches the wall, the longitudinal DOS function in- 
creases in a similar way as the transverse DOS function 
does in Fig. 0(a). The vertical dashed lines indicate the 
'bare' atomic transition frequencies for which the atomic 
vdW energies in Fig. 0(b) were calculated. Figure 0(b) 
represents the exact solutions of Eq. (|48|l for atom inside 
and outside the (10,10) CN. When the atom is inside, the 
vdW energies are seen to be in general lower than those 
for the atom outside the CN. This may be attributed to 
the fact that, due to the nanotube curvature, the effective 
interaction area between the atom and the CN surface is 
larger when the atom is inside rather than when it is out- 
side the CN. This, in turn, indicates that encapsulation 
of doped atoms into the nanotube is energetically more 
favorable than their outside adsorption by the nanotube 
surface — the effect observed experimentally in Ref. |5j. 

Comparing Fig. (b) with Fig. (b) for the outside 
atomic position, one can see that, if the atom-CN-surface 
distance is so large that the weak-coupling approximation 
is good, then for the same distance and approximately 
the same atomic transition frequency xa ~ 0.3 the vdW 
energy near the (10,10) nanotube (R cn = 6.78 A) is lower 
than that near the (9,0) nanotube (R cn — 3.52 A). This 
is in agreement with the qualitative predictions formu- 
lated above in analyzing the weak coupling regime. In 
the present case, the effective atom-CN-surface interac- 
tion area is larger for the atom outside the larger-radius 
(10,10) nanotube, thereby explaining the effect. Obvi- 
ously, one should expect an opposite effect for the inside 
atomic position. Experimentally, the heat of absorbtion 
of Krypton, to take an example, has been shown to be 
larger for graphite (0.17 eV/atom) than for multiwall 
nanotubes (0.12 eV/atom) 30], in agreement with our 
theoretical predictions. 

A similar tendency is demonstrated in Fig0 Here, the 
vdW energies calculated from Eqs. (^HJ-JSUJ are repre- 
sented for the atom positioned at a fixed distance from 
the nanotube wall [chosen to be equal to the radius of 
the (9,0) CN] inside and outside the "zigzag" (m, 0) CNs 
of increasing radius R cn — m (ay3/2-K) (a = 1.42 A is the 
CN interatomic distance). When the atom is inside, the 
energy absolute value goes down with R cni representing 
a general tendency of the effective interaction area reduc- 
tion with lowering the CN surface curvature. When the 
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FIG. 5: Ground-state van der Waals energy of the two-level 
atom positioned at a fixed atom-surface distance [equal to the 
radius of the (9,0) CN] inside and outside "zigzag" (m, 0) CNs 
of different radii, as a function of the nanotube index m. 



atom is outside, the effect would be the opposite one if 
the atom were weakly coupled to the field, thus indicat- 
ing that the atom-field coupling regime is not actually 
weak for the distance chosen. In general, the inside and 
outside vdW energies approach each other with R cn as it 
should be since the two atomic positions become equiv- 
alent in the plane limit where R cn — > oo. Interesting is 
also the fact that the vdW energies of the atom outside 
the metallic nanotubes (m is divisible by 3) are on aver- 
age a little bit smaller than those for the atom outside 
semiconducting ones. This is certainly the property re- 
lated to the difference in the conductivities of metallic 
and semiconducting "zigzag" nanotubes. 



V. CONCLUSIONS 

Recently, a theory of atomic spontaneous decay near 
CNs was developed 0,0! El - The theory has brought out 
fascinating peculiarities of the vacuum-field interactions 
in atomically doped CNs. In particular [a,E|, *f the atom 
is close enough to the CN surface and the atomic tran- 
sition frequency is in the vicinity of a resonance of the 
photonic DOS, the atomic spontaneous decay dynamics 
exhibits vacuum-field Rabi oscillations. This is a prin- 
cipal signature of strong coupling of the excited atomic 
state to vacuum surface photonic modes in the CN. Such 
a strongly coupled atomic state is nothing but a 'quasi- 
1D cavity polariton' similar to quasi-OD excitonic (elec- 
tronic) polaritons in quantum dot microcavities widely 
discussed in the literature |3li |32j, |33J . The study of such 
phenomena was started awhile ago in atomic physics |34j 
and still attracts a lot of interest especially in connection 
with research on quantum computation [35| . 

Obviously, the stability of quasi-lD atomic polaritons 
in CNs is mainly determined by the atom-nanotube van 
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der Waals interaction. This interaction, however, orig- 
inates from strong atom-field coupling and, therefore, 
cannot be correctly described in terms of vacuum-QED- 
based (weak-coupling) vdW interaction models as well 
as in terms of those based upon the linear response the- 
ory. In the present paper, we have developed a simple 
quantum mechanical approach to the ground-state vdW 
energy calculation of a (two-level) atomic system near a 
CN. The approach is based upon the perturbation the- 
ory for degenerated atomic levels, thus accounting for 
both weak and strong atom-field coupling and thereby 
covering the vdW interactions of the ground-state quasi- 
1D atomic polaritons in CNs. Within the framework of 
this approach, the vdW energy is described by the in- 
tegral equation represented in terms of the local pho- 
tonic DOS. By solving this equation numerically, we have 
demonstrated the inapplicability of weak-coupling-based 
vdW interaction models in a close vicinity of the nan- 
otube surface where the local photonic DOS effectively 
increases, giving rise to an atom-field coupling enhance- 
ment followed by the degeneracy of the unperturbed 
atomic levels due to the diamagnetic interaction term. 

The fundamental conclusion above is supplemented 
by those important for various applications of atomi- 
cally doped CNs in modern nanotechnology. In partic- 
ular, we have studied CN surface curvature effects on the 
atom-nanotube vdW interactions and have shown that 
an inside encapsulation of doped atoms into the nan- 
otube is energetically more favorable than their outside 
adsorption by the nanotube surface, in agreement with 
the experimental observations reported in Ref. [f| . More- 
over, if the atom is fixed outside the CN in a way that 
the weak atom-field coupling regime is realized, i.e. far 
enough from the CN surface, then the modulus of the 
atom-nanotube vdW energy increases with the CN radius 
because the effective atom-nanotube interaction area is 
larger for larger-radius nanotubes in this case. For inside 
atomic position, the modulus of the atom-nanotube vdW 
energy decreases with the CN radius, representing a gen- 
eral effect of the effective interaction area reduction with 
lowering the CN surface curvature. 

A further intriguing extension of the present work 
could be the vdW interactions of excited atomic states 
where, even in the weak atom-field-coupling regime and 
in the simplest case of an atom near a planar semi-infinite 
medium, very interesting peculiarities (e.g., oscillatory 
behavior) were recently shown to exist [12 ] . 
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APPENDIX A: THE TWO-LEVEL 
APPROXIMATION 

In this Appendix, we rewrite the Hamiltonian i|19|) - 
((23(1 in terms of a two-level atomic model [24j. Within 
the framework of this model, the spectrum of the atomic 
Hamiltonian ((21(1 is approximated by the two eigenstates, 
upper \u) and lower \l), with the energies Hlo u and htoi, 
respectively. One has then 



H A w hw u \u)(u\ + fiwi\l){l\ 
with a completeness relation 

|u)(u| + |/)<i| = J. 



(Al) 
(A2) 



Subtracting a constant term (h/2)(u) u + cj;)7 from the 
Hamiltonian IjAlfl and thereby placing the energy zero in 
the middle between the two energy levels, one arrives at 
the 'bare' two-level atomic Hamiltonian of the form 



H 



hujA 



(A3) 



with a z = \u) (u\ — \l) (l\ and lo a = uj u — uji being the 'bare' 
atomic transition frequency. 

In terms of such a two-level scheme, the atomic dipole 
moment operator d has the matrix elements (u|d|Z) — 
(l\d\u) — d and (u|d|tt) = (Zj_d.|Z) = 0. Keeping this and 
the completeness relation ((A2(l in mind, the interaction 
Hamiltonians (12211 and 1(23(1 can be rewritten as follows. 

For the interaction J22J, using Eqs. © and (fTTTft and 
a well-known quantum mechanical operator equality 

(pi) a = mi—(ri) a = ^ [(Ti) a ,H A ] , (A4) 
at ih 

where the Latin and Greek indexes enumerate particles 
in the atom and vector components, respectively, and Ha 
is given by Eq. ((A3() . one has 



AF 



(a + a^d- 



aw — E (ta, — h.c. 



cLuE (ta, uj) + h.c. 



(A5) 



with a — \l)(u\ and Sfi — Using further Eqs. ltTT|l . 

()17J) and (0) and assuming the longitudinal (along the CN 
axis) atomic dipole orientation due to the dominant axial 
polarizability of the nanotube [l^, Hil j one arrives 

at the secondly quantized interaction Hamiltonian l|27() 
with the interaction matrix elements 130(1 and l(31|l . 
For the interaction 1(23(1 . one may proceed as follows 



A a (r A )Ap(r A ) 



2hc 2 



y2^[(Pi) a ,dp}A a (rA)A (r A ), (A6) 
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where the product of the two Kronecker-symbols in the and the completeness relation (|A2|) . one arrives at 
first line is represented in terms of the 'coordinate- 



momentum' commutator as 



H 



(2) 
AF 



Si 3 S a (3 = - [(Pi) a , (Tj)p]. 



UJA_ 

he 2 



d-A(r A ) 



(A7) 



which, upon substituting Eq. Q for the vector potential 
Using Eq. ((A4() for the momentum operator components operator, becomes 



du> 



LUA -L -L f ^ A f ^ A -L r 



h.c. 



r 



The four items here can, in view of Eqs. (|ll(l . 117(1 and 

(0, be classified as follows. The first item and its her- 
mitian conjugate describe the processes with simultane- 
ous annihilation and creation of two photons (two-photon 
transitions). They may be safely neglected if the field in- 
tensity is low enough that, we believe, is the case for the 
vacuum electromagnetic field we deal with in consider- 
ing the vdW interactions. Furthermore, keeping them in 



~ (2) 

the Hamiltonian H A p is meaningless as the Hamiltonian 

H AF is nothing but a dipole approximation which, by its 
definition, neglects two-photon transitions. The second 
item and its hermitian conjugate are combined together 
by using bosonic commutation relations (J2J to give in 
terms of the bosonic field operators and the transverse 
dipole interaction matrix element l(31|) the expression 



p- /l/dR|g^(r A ,R,c)| 2 -|^ 
riuj A J J Tvuj A 



flu du'fdRdR' [g^, R, w)] V(r A , R', ^^(R w)/(R', u>% 



r 



where the second term is nothing but a two-photon cor- 
rection to the dipole interaction 1(27|). that can be easily 
seen by noting that a z = a^a — aa^ . This correction must 
be neglected for the reasons just discussed, so that finally 
one arrives at 



H 



(2) 
AF 



HuJ A 



da; dR\g^(r A ,R,Lj)y 



(A8) 



Eq. (|A8|I can now be combined with the 'bare' atomic 
Hamiltonian i(A3(l to yield the 'effective' unperturbed 
atomic Hamiltonian 1)26(1 with the renormalized atomic 
transition frequency given by Eq. ((29(1 . 

Thus, we have proved that the Hamiltonian ((19(l - ((23(l . 
being approximated in terms of the two-level atomic 
model, is represented by the Hamiltonian (|21|) - (|T7I) . 



APPENDIX B: THE TRANSVERSE AND 
LONGITUDINAL LOCAL PHOTONIC DOS 

We start with Eq. 1(18(1 for the Green tensor of the elec- 
tromagnetic subsystem. This equation is a direct conse- 
quence of the equation 

A-7T 

£ (VxVx-fc 2 ) 2a | a (r,,)^-4(r, W ) (Bl) 



obtained by substituting the magnetic field operator from 
Eq. lJT3jl into Eq. (|T5|l . On the other hand, under the 
condition of the Coulomb gauge (see, e.g., |3(|), one has 

E(r,o;) =E" L (r,o;) +E (r,w) (B2) 
with [compare with Eqs. © and tfTT3f> ] 

E X (r,w) =ikA(r,u), (B3) 

where 

V-A(r,cj) = 0, (B4) 

and 

E"(r,w) = -Vip(r,(j). (B5) 

From Eqs. ((B1(| - [(B5() . in view of Eq. {llj and {T7|), one 

obtains the following equation for the Green tensor com- 
ponents 

^(VxVx- fc 2 ) J x G Q ,(r, R, u) + ^G az (r, R, u 



a—r,ip,z 



= a(r-R) 
with additional constraints 

V Q x G QZ (r,R,^) = 



(B6) 
(B7) 



a— r,(p,z 
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and 

ea/3 7 V,gllG 7Z (r,R,cj) = 0, (B8) 

where e Q/ 3 7 is the totally antisymmetric unit tensor of 
rank 3. Keeping Eqs. ifBTJl and (|B8)l in mind, Eq. l|B6|) 
is rewritten to give two independent equations for the 
transverse and longitudinal Green tensors of the form 

(A + fc 2 ) x G zz (r, R, w) = -6i(r - R) , (B9) 
fc 2 llG zz (r,R, W ) =-4(r-R) (BIO) 

with the delta-functions given by Eqs. (|12|) and i|13|) and 
±( -^G zz defined by Eq. (^J. 

The transverse Green function ± G zz (r, rj, oj) was de- 
rived and analyzed in our previous work [lOj . Its differ- 
ential form equivalent to Eq. i|32[l is given by 

± G 2Z (r, r A , lj) = (^V Z V Z + 1 j g(r, r A ,u) (Bll) 

with g(r, r A ,u>) being the Green function of the scalar 
Helmholtz equation, satisfying the radiation condition at 



infinity and boundary conditions on the CN surface. In- 
deed, substituting Eq. IBllfl into Eq. I|B9[) and using 
Eq. (|BTfl , one straightforwardly obtains 

(A + fc 2 ) g(r, r A , to) = -<^(r - v A ) , (B12) 

the scalar Helmholtz equation with a transverse ^-source. 
For the longitudinal Green function "G zz (r, Ta, w), one 
analogously has 

"G„(r,rA,w) = —j3V z V x g(T,T A ,u), (B13) 
which upon substituting into Eq. l|B10fl yields 

V z V z g(r,r A ,u) = 6l(r-r A ). (B14) 

The complete solution of Eq. (|B11|) , which satisfies the 
radiation condition at infinity and boundary conditions 
on the CN surface, has for r A > r > R cn the form 0] 



g{r,r Ai u) =g {r,r A ,u)) - — ^ > e w / ^ 

( 27F ) pf^oo "/C 1 + £H V 2 RcnIp{V, 



p (wr J 4)iv: p (wr) 

Jp(vR cn )Kp(vR cn ) 



e ihz dh , 



(B15) 



where the coordinate system has been fixed as is shown 
in Fig. ^ so that r A = {r A ,0, 0} and r = {r, <p,z}, the 
first item is the particular solution of the inhomogeneous 
equation which is 



#o(r,r A ,w) 



^ gifc|r-rA 

4tt |r - r A \ 



(B16) 



the point-source function (see, e.g., Ref. 14]), and the 
second item is the general solution of the homogeneous 
equation with the boundary conditions imposed on the 
CN surface and at infinity. This latter item is repre- 
sented in terms of the modified cylindric Bcssel func- 
tions I p and K p with v = v(h,oj) = \/h 2 — fc 2 and 
f3(u>) — Airi(j zz (R cn , ui)/u>. The integration contour C 
goes along the real axis of the complex plane and en- 
velopes the branch points ±fc of the integrand from below 



J 



and from above, respectively. The function g(r, r A , ui) for 
r < r A < R cn is obtained from Eq. IJB15II by means of a 
simple symbol replacement I p «-> K p in the numerator of 
the integrand. 

Now, using Eqs. (|B1 1|) . I|B13|I . (|B7|) and the property 
g(r, r A , ui) = g(r A , r, uj) which is obvious from Eqs. (|B15|) 
and (|B 16|) . it is not difficult to prove that 

± Gj z (r A ,r A ,u;) = x G zz (r A , t a ,u) , (B17) 
^GI z (ta,t a ,lj) = ^G tz (T Ai T A ,w), (B18) 

(B19) 



X GI{T A ,T A ,U) = 0, 



whereupon, making use of the definitions l|35|) . i|36[l and 
Eqs. I|38|) - (|42|l . one arrives for r A > R cn at the the trans- 
verse and longitudinal local photonic DOS functions of 
the form 



oo 



7:-L / \ 3i? C n T 



dh[3{u)v 4 I 2 (vR cn )K 2 (vr A ) 



27rfc 3 p £^ o J c 1 + I3(lu) v 2 R cn Ip(vR cn )K p {vR cn ) ' 

2^ Im ^ J c l 



p— — oo 



c 1 + /?M v 2 R cn Ip{vR cn )K p (vR cn ) 



(B20) 
(B21) 
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Here, the sign in front of Eq. I(B21(I has been changed to 
be positive. This reflects the fact that the right-hand sides 
of Eqs. (|B12|) and (|B14(I are of opposite signs. As a con- 
sequence, partial solutions of the corresponding homoge- 
neous equations should be taken to have opposite signs 
as well. This yields a correct (positive) sign of the longi- 
tudinal local photonic DOS which, along with the trans- 
verse local photonic DOS, must be a positively defined 
function. For ta < R cn , Eqs. (|B20|) and (|B21|) should be 
modified by the replacement ta *-* R C n m the Besscl func- 
tion arguments in the numerators of the integrands. In 
dimensionless variables l|47|) . these equations are rewrit- 
ten to give Eqs. lO and (EOJ. 

Finally, based upon the properties (|B 1 7(1 - (|B 19(1 . one 
can easily prove Eqs. ((33(1 and 1(34(1 for the dipole inter- 
action matrix elements. The proof of Eq. 1(33(1 is straight- 
forward. In view of Eq. 1(31(1 . one has 

/I R/j2 J 2 r 
dR\g ± ^(r A ,R,Lu)\ 2 =irhcu— f-i J dRRea zz (R,u) 

x J -(lliG M (r A ,R,w)- L (IDG„(r A ,R,w)*. (B22) 
Using further the integral relationship 

47T 

ImG Q/3 (r, r',Lu) = — k 
c 

x J dRRea zz (R,u;)G az (r,R,Gj)Gf} z (r',R,(jj (B23) 

[which is nothing but a particular 2D case of a general 
integral relationship proven for any 3D electromagnetic 
field Green tensor in Ref. 0], with Eq. taken into 
account] , one obtains for the right-hand side of Eq. ((B22() 

n ^l ln ^w) G n\\ ){YAiVA ^ 

c 

whereupon using Eqs. (j3"5*|l . ((37))) and (jSHJ, one arrives at 
Eq. 123). 



The proof of Eq. ((34(1 is based upon Eq. 133(1 . In view 
of this equation, the right-hand side of Eq. ((34(1 takes the 
form 



JdR\g^(r A ,R,cu)\ 2 



(B24) 



± — 2Re / dR g 1 - (r A , R, lo) gU (r A , R, wf, 
u A J 



where the last item can be further rewritten in terms of 
Eqs. (EIJ), (E2J> and l|B23)l to give 



which is zero, according to Eq. ((B19(l . Thus, one arrives 
at Eq. PHl. 



APPENDIX C: THE RELATION OF EQ.lHSl) 
WITH THE CASIMIR-POLDER FORMULA 



In deriving an 'infinitely conducting plane' result (the 
Casimir-Polder formula 27]) from our theory, we will fol- 
low a general line of the work by Marvin and Tiogo [37j 
who did the same within the framework of their linear 
response theory. 

We start with our weak-coupling-regime equation 1(54(1 . 
First of all one has to simplify the local DOS functions 

£ (r A ,x) in this equation by putting ~a zz — ► oo. From 
Eqs. ((49(1 and ((50|) . one has then for ta > R C n 



£ (r A ,x) 
I (r A ,x) 



— Im 



dy 



y — 1 — ie 

y 2 



E 



KplVv 2 - 1 i£ u(r A )x}I p [y/y 2 - 1 - ie u(R cn )x] 
s K p [x^Jy 2 - 1 - ieu{R cn )x] 



r 



(Cl) 



where e is an infinitesimal positive constant which is nec- 
essary to correctly envelope the branch points of the in- 
tegrand in integrating over y. 

The next step is taking a large radius limit in Eq. I(C1(I . 
This can be done by using the relationship 



lira 

a — > oo n - 
b — > oo 
(a — b — const) 



E 



K 2 (a)Ub) 
K n {b) 



K [2(a-b)] 



proved in Ref. |3j| ■ One has 
fl ± (r A ,x)\_ 3 



/•oo 

Im / dy 

Ti" Jo 



1 — ie 



x K [2fix\/y 2 -l-ie 



(C2) 



where fi = 2"fol/hc with I = r A — R C n being the atom- 
surface distance. The integration variable y in this equa- 
tion contains implicit frequency dependence. Indeed, 
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comparing dimensionless equations (|49(l , (|50(l with (|B20|) , 
(|lj21|) . one can see that y — h/k — h(c/uj) = h(ch/2j x). 
For the following it makes sense to extract this frequency 
dependence by the substitution y — r/x with dimension- 
less t — h(cTi/2^o). The result takes the form 



with 



T M\\) 



(ta,x) = ^ {f^(r Al x) -f^(r A ,xf} 



(C3) 



x K [2px\J (t/x) 2 — 1 — is ] . 
Substituting Eq. (|C3(1 into Eq. l(5i|l. one has 

3x A 



(C4) 



4l7T 2 



<ix Tq (a;) 



(lA+l)l 2 

dx To (a;) 



f- L (r A ,x)-—p(r A ,x) 
x A 



(C5) 



o (au-ai)x 2 



/ x (r A ,-a;) +— /"(r Al -z) 



Here, in the second item, the change of the integration 
variable from x to — x has been made and it has been 
taken into account that Tq(—x) — — Fq(x). Noting fur- 
ther that both integrands in Eq. (|C5|) do not have poles 
in the upper complex half plane, one may rotate the in- 
tegration path of the first and the second integral by ir/2 
and by — 7r/2, respectively. Then, both integrals become 
the integrals over the positive imaginary axis of the com- 
plex plane and Eq. I|C5(I takes the form 



3x A 

4l7T 2 



dz Tq(z) 



o (x A + z)z 2 
dz Tq(z) 



f X (r A ,z) fHr A ,z) 

x A 



(C6) 



(x A - z)z 1 



f ± (r A ,-zf + — fHr A ,-zf 
x A 



Making further the change of the integration variable 
from z to iu with real non-negative u, one arrives at 



£vw(ya) 



3x A 
An 2 




du To(iu) 



o (x A +iu)u 2 
duTo(iu) 



f^{r A ,iu) p(r A ,iu) 

x A 



(C7) 



(Xj 



m f\\t ■ f 

P(r A ,-m) 

x A 



Here, the functions / _L ^^(r J 4, iu) have the following ex- 
plicit form 



'f ± {rA, iu) 
ftt(r A ,iu) 



h I VVT - 1 



if (2 MW? )(C8) 



obtained from Eq. 1|C4|I after the change of the integra- 
tion variable by the subsequent substitutions x = iu and 
?7 2 =t 2 /u 2 + 1. The functions f^\-iu)*=f^\iu)as 
a consequence of a general property K v {z)* —K v {z*) [26j. 
Substituting Eq. l|C8f> into Eq. HC7f) and making al- 
lowance for the fact that Tq(iu) — — iTo(u), one further 
has 



du f o(u) 




which upon the substitution r/ 2 — p becomes 



K (2uur)) , (C9) 



/ \ 3x A 



' du To (it) 
x 2 A + u 2 



^dpL/p~^l + -^=^K (2au^p). (CIO) 



In Eq. (jCTU|) . the first integral over p is taken exactly 
by means of the relationship 38] 



dp\J p- lK (ay/p) = \j if 3/2 (a) 



:(l + a)e- a 



and the second one - by making a large argument series 
expansion of the integrand and its subsequent termwise 
integration. This, upon returning back to dimensional 
variables by means of Eqs. I|47fl and using the definition 
(see, e.g., HJ) 



a zz (iu) 



d z uj A 
\-{iuf 



for the polarizability tensor zz-component of the two- 
level system, finally yields 



E vw (r A ) 



8ttP Jo 



dua zz (iu) 



-2{u/c)l 



This is the half of the well-known Casimir-Polder result 
for the vdW energy of an atom near an infinitely con- 
ducting plane (see Refs. H3,|33|), in agreement with the 
fact that our model of the atom-electromagnetic-field in- 
teractions in the presence of a nanotube only takes the 
longitudinal (along the nanotube axis, or, equivalently, 
parallel to the plane in the large nanotube radius limit) 
atomic dipole orientation into account. Another half of 
the Casimir-Polder vdW energy would come from the 
transverse (perpendicular to the plane) atomic dipole ori- 
entation which we have neglected. 
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